Realization of Berezinskii's superconductivity in quasi-one-dimensional systems 
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We revisit the pairing symmetry competition in quasi-one-dimensional systems. We show that 
spin-triplet s-wave pairing, where the pair is formed by electrons with different times and has an odd- 
frequency symmetry, can be realized in systems with strong one-dimensionality when the strength of 
charge fluctuation dominates over spin fluctuation. The present study provides a novel microscopic 
mechanism for this exotic pairing originally proposed by Berezinskii in 1974. 



Superconductivity in strongly correlated systems has 
been a long standing issue in condensed matter physics. 
To avoid the strong on-site Coulomb repulsion, a Cooper 
pair tends to be formed by two electrons located in sep- 
arate places with the non-zero angular momentum L. 
Spin-singlet d-wave (L = 2) [1| and spin-triplet p-wave 
(L = 1) [2| pairings belong to this class. On the other 
hand, two electrons located in separate times are also 
able to avoid the Coulomb repulsion in forming a Cooper 
pair. There is a possibility that a sign of the pairing 
function changes in exchanging times of two electrons. 
Then the Fourier transformed pairing function is an odd 
function of frequency. This type of pairing is called odd- 
frequency (odd-w) one originally proposed by Berezinskii 
[3|, while conventional pairing is called even- frequency 
(even-w) one. 

Berezinskii proposed odd-w spin-triplet s-wave pairing 
(Berezinskii's pairing) in a bulk system. However, real- 
ization of such a pairing has been believed to be difficult 
in bulk systems due to the presence of gapless excitation. 
Bergeret et al. revived this exotic pairing in the context 
of ferromagnet / superconductor (F/S) issues |j|. They 
have proposed that Berezinskii's pairing can be induced 
locally in ferromagnet although the bulk pairing symme- 
try in superconductor is spin-singlet s-wave. Stimulated 
by this proposal there have been several works about this 
pairing in F/S junctions [5ll6|. Besides ferromagnet junc- 
tions, Berezinskii's pairing has been predicated in diffu- 
sive normal metal attached to spin-triplet p-wave super- 
conductor [3]. It is a really interesting issue to resolve 
whether Berezinskii's pairing is possible in bulk super- 
conductor or not. 

Strictly speaking, there are two classes of odd-w super- 
conductors, i.e., spin-singlet and spin-triplet ones. Odd- 
Lo spin-singlet superconductor _&] has an odd-parity in ac- 
cordance with Fermi-Dirac statistics. Thus the resulting 
superconducting state is fragile against impurity scatter- 
ing similar to spin-singlet d-wave and spin-triplet p-wave 
pairing ones. Furthermore, temperature dependence of 
Knight shift below transition temperature Tq becomes 
similar to that of spin-singlet d-wave superconductor. On 
the other hand, odd-oj spin-triplet s-wave superconduc- 
tor, the original proposal by Berezinskii, can have a clear 



difference from preexisting pairings, since it has an un- 
changed Knight shift below Tq even in the presence of 
impurity scattering. Thus, in order to discover odd-w 
superconductors experimentally, study on the mechanism 
of Berezinskii's pairing is highly desired. Although there 
have been several studies on the generation of odd-w 
pairing up to now |8l-ll0l|. there has not been clear mi- 
croscopic mechanism which supports the realization of 
Berezinskii's pairing. The aim of the present paper is to 
present a clear microscopic mechanism for realizing this 
exotic pairing. 

To resolve this issue, we focus on quasi-one- 
dimensional (QID) systems. Up to now, superconduc- 
tivity in QID systems has been studied in the context 
of, e.g., an organic superconductor (TMTSF)2X. It has 
been shown that spin-triplet /-wave pairing (L = 3) can 
dominate over spin-singlet d-wave {L = 2) one [ll| when 
charge fluctuation dominates over spin fluctuation |12l - 



14| as we shall discuss below. In real space, these corre- 



spond to Cooper pairs formed between separate places. 
However, when the system becomes strongly QID, it is 
difficult to form Cooper pairs with separate places due to 
the geometrical constraint, so actually there is a chance 
that even- and odd-w pairings compete. In the present 
paper, combining this "odd-w > even-oj" situation with 
the "triplet > singlet" effect, we show that spin-triplet 
s-wave {L — 0) pairing, namely, the original Berezinskii's 
pairing, can be realized in a strongly QID system with 
strong charge fluctuation. This is actually exemplified by 
solving the linearized Eliashberg's equation in the QID 
extended Hubbard model. The present study provides 
a novel and realistic mechanism for realizing this exotic 
pairing proposed more than thirty years ago. 

Before going into the actual model and the calcula- 
tion results, we make a general argument for the pairing 
symmetry in a QID system. We assume a many body 
system on a QID lattice, where the hopping integral in 
the y direction ty is smaller than that in the x direction 
tx. The on-site Coulomb repulsion enhances spin fluctu- 
ation at the nesting vector Q. In addition, we assume a 
situation where the off-site Coulomb repulsion enhances 
charge fluctuation at Q. When the pairing interaction 
is mainly mediated by spin and charge fluctuations at 



Q, the effective pairing interactions for spin-singlet and 
spin-triplet channels can be given by 
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Kff(il^m,Q) = -2^sp(ii^m,Q) - 2^ch(ii^m,Q), (2) 

respectively, where T^p and Vch are contributions from 
spin and charge fluctuations, respectively. Vm — 2rmrT 
is the bosonic Matsubara frequency with an integer m 
at the temperature T. In strongly correlated systems, 
the effective pairing interaction at i/m = tends to give a 
large contribution to pairing. When the off-site Coulomb 
repulsion is absent or small, spin-singlet pairing is favored 
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due to the prominence of spin fluctua- 
tion. We call this "SF > CF case" . On the other hand, 
when the off-site Coulomb repulsion is so remarkable as 
to make charge fluctuation exceed spin fluctuation, spin- 



triplet pairing is favored (|V^ff| > 
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"CF > SF case" . It is also noted that the effective pair- 
ing interaction for spin-singlet (spin-triplet) channel has 
a positive (negative) sign as far as the contribution by 
charge fluctuation does not become too large. Since the 
superconducting gap function A has to satisfy a condi- 
tion ^•\Q)A{kF)A{kF + Q) < on the Fermi surface, it 
is required for spin-singlet (spin-triplet) pairing to satisfy 
a condition A(fcF)A(fcF + Q) < (A(A:f)A(A:f + Q) > 0) 
with fcF — (ia;„,fcF) and Q = (0,Q) consisting of the 
momentum on the Fermi surface fcF and the fcrmionic 
Matsubara frequency w„ = (2?! — IjirT. In the follow- 
ing, we discuss all four classes of pairings in accordance 
with the Fermi-Dirac statistics. They are (i) even-w spin- 
singlet even-parity (ESE), (ii) even-w spin-triplet odd- 
parity (ETO), (iii) odd-w spin-singlet odd-parity (OSO), 
and (iv) odd-w spin-triplet even-parity (OTE) pairings. 
Berezinskii's pairing belongs to class (iv). We consider 
four cases by combining strongly (weakly) QID lattice 
and SF > CF (CF > SF) case. 

First, we discuss a weakly QID lattice. In this case, 
electrons avoid the strong Coulomb repulsion in real 
space. Anisotropic pairing with the non-zero angular 
momentum is induced and it has even-w symmetry. It 
is difficult for odd-w pairing to be stabilized due to the 
nature of pairing with separate times. In the SF > CF 
case, which favors spin-singlet pairing, ESE one has an 
advantage. This pairing has d-wave symmetry, where two 
nodes of A run close to the Fermi surface as shown in Fig. 
[T](A). Reversely in the CF > SF case, which favors spin- 
triplet pairing, ETO one dominates. This pairing has 
/-wave symmetry, where two nodes of A also run close 
to the Fermi surface as shown in Fig. [l](B) |15j . 

Next, we consider a strongly QID lattice. In this case, 
avoidance of electrons is limited in real space. This can 
induce pairing with separate times and make odd-w pair- 
ing comparable to even-w one. If we focus on the above 
leading pairings d- and /-wave, the gap nodes will run 




FIG. 1: (color online). Pictures of the momentum k = 
(ka:,ky) dependences of A. Upper panels show sign change 
(positive or negative) of A with nodes (black dashed lines) 
in momentum space. (Purple) solid curves denote the Fermi 
surface with the nesting vectors (arrows). Lower panels show 
the kx dependences of A. (A) d-wave in SF > CF case on a 
weakly QID lattice. (B) /-wave in CF > SF case on a weakly 
QID lattice. (C) p-wave in SF > CF case on a strongly QID 
lattice. (D) s-wave in CF > SF case on a strongly QID lattice. 



almost on the entire Fermi surface (see Fig. [T](A),(B)), 
which destabilizes those pairings. Odd-o; symmetry al- 
lows pairing with lower angular momentum than that 
of even-w one, namely, the number of gap nodes in mo- 
mentum space is smaller. As a result of the competition 
between gap nodes in momentum and frequency spaces 
(see the insets of Fig. ^ , odd-w pairing can replace even- 
uj one as leading one. In the SF > CF case, OSO pairing 
has an advantage. This pairing has p-wave symmetry, 
which has no nodes of A on the Fermi surface in momen- 
tum space as shown in Fig. [1] (C). Finally in the CF > 
SF case, spin-triplet pairing dominates over spin-singlet 
one, and the OTE pairing can take place. This is in- 
deed the spin-triplet s-wave pairing originally proposed 
by Berezinskii. Here, note that the mechanism is en- 
tirely novel in that the odd-w > even-o; situation due to 
the strong one-dimensionality combined with the triplet 
> singlet situation given by the CF > SF case is the ori- 
gin of the realization of this exotic fully gapped state in 
momentum space (Fig. [I](D)). 

In order to exemplify the above physics, we actually 
apply the random phase approximation (RPA) to the 
QID extended Hubbard model and solve the linearized 
Eliashberg's equation. The Hamiltonian considered here 
is given by 



n = -Y.t 
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tij is the hopping integral between sites i and j. The 
hopping integrals between sites neighboring in the x, y, 
and diagonal directions are labeled as i^, ty, and id, re- 
spectively as shown in Fig. [H We take account of the 
diagonal hopping integral that enhances the geometri- 
cal frustration, which can significantly induce competi- 
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FIG. 2: (color online). The model considered in the present 
study, tx, ty, and td are the hopping integrals. U and V are 
the on-site and off-site Coulomb repulsions, respectively. 



tion between even- and odd-w pairings. The dispersion 
is given by Efc — ~'2txCoskx~2tyCosky — 2tdCOs{kx + ky). 
We choose ty =td — O.SSi^: {ty —td — O.lij^) as the hop- 
ping integrals for a weakly (strongly) QID lattice. c^J 
and riic are the annihilation (creation) and number oper- 
ators for an electron with spin cr on a site i. rii = rii-^+nii. 
U and V are the on-site and off-site Coulomb repulsions, 
respectively. V acts between electrons neighboring in 
the X direction. The momentum dependence is given by 
V{q) — 2T^cosga; with the momentum q = {qx^Qy)- In 
this model, we solve the linearized Eliashberg's equation 
for spin-singlet (spin-triplet) channel within the RPA 



^^(^) = -^ E ^cf (fc - fc')|Go(fc')pA(fc'), 

k' 



(4) 



where N is the number of sites and k = (ia;„, fc). G'o(/c) = 
(iujn — Efe + m)""'^ is the bare Green's function with the 
chemical potential fi. A is the eigenvalue for A. A be- 
comes unity just at Tc- The more stable the supercon- 
ducting state is, the larger A tends to be. We calcu- 
late A and A for ESE, ETO, OSO, and OTE symme- 
tries. The effective pairing interactions for spin-singlet 
and spin-triplet channels within the RPA are given by 

,{q) = U + V{q) + ^U\sp{q) 

-l{U + 2V{q)}\,i,{q) (5) 
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respectively, with q = {ivm,q)- 
susceptibilities are given by 

xo(g) 



2^^(g)}'xch(g), (6) 

The spin and charge 



Xsp(g) = 
Xch(g) == 



1 



Uxoiq) 

xo{q) 



l + {U + 2V{q)}xo{q)' 
respectively, with the irreducible susceptibility 
T v-^ 



xoiq) 



k 



Go{q + k)Go{k). 



(7) 
(8) 

(9) 



In the present paper, we normalize the gap function 
J2k 1^(^)1/-^ = 1- We take the number of sites N = 
Nx X Ny = 256 X 64. The number of electrons is 
unity per site (half- filled) , which gives the nesting vec- 
tor Q — (7r,7r/2) both on weakly and strongly QID lat- 
tices. The fermionic and bosonic Matsubara frequency 
have values from -(2A'^f - l)7rT to (2A'^f - l)7rr and from 
-2NiTTT to 2NiTTT, respectively, with Nf = 2048. 
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FIG. 3: (color online). T dependences of A in SF > CF (CF 
> SF) case on a weakly (strongly) QID lattice. Insets (k) and 
(w) are fc = {kx, ky) and u)n dependences of A corresponding 
to the largest eigenvalue, respectively. In the inset (k), we 
choose ky = G and ujn = ttT. /cf in the inset denotes the 
location of the Fermi surface at ky = 0. In the inset {u)), we 



choose fe where A gets the largest value on ky 



0. 



In Fig. [3l A for four kinds of pairings are plotted 
against T in the four cases. There are two insets in pan- 
els Fig. [31JA)-(D), where (k) k^ and (w) w„ dependences 
of A are plotted at T = O.Obt^. In the inset (k), we 
choose ky = and w„ = nT. k-p in the inset denotes 
the location of the Fermi surface at ky — 0. In order to 
give clear comparison with inset (w), we plot only the 
fca; > portion, while the fc^^ < portion is given by 
A(A;) = +(— )A(— fc) for even- (odd-) parity pairings. In 
the inset (w), we choose k where A gets the largest value 
on ky = 0. 

First, we focus on a weakly QID lattice [ty — td ~ 
0.35U). In the SF > CF case {U = 2^, V = 0), ESE 
pairing is the most stable as shown in Fig. EIJA). This 
pairing has d-wave gap, whose fc^: > portion is shown 
in the inset (k) of Fig. [31[A). There is a node of A near 
kx = kp in momentum space, while there are no nodes in 
Matsubara frequency space as shown in the inset (w) of 
Fig. M.A)- In the CF > SF case {U = 1.995^^^, V = t^), 
ETO pairing is the most stable as shown in Fig. El^B). 
This pairing has /-wave gap, whose kx > portion is 
shown in the inset (k) of Fig. |3l^B) [12|. There is also a 
node of A near kx = k-p in momentum space, while there 
are no nodes in Matsubara frequency space as shown in 
the inset (w) of Fig. EjB). 

Next, we focus on a strongly QID lattice [ty ~ td ~ 
O.ltx). In the SF > CF case ([/ = 1.6t^, V = 0), OSO 
pairing is the most stable as shown in Fig. [3l^C). This 
pairing has p-wave gap, whose fc^; > portion is shown 
in the inset (k) of Fig. [3lJC) [10|. There are no nodes of 
A near kx = kp in momentum space, and instead there 
is a node in Matsubara frequency space as shown in the 
inset (w) of Fig. ^C). And finally for the CF > SF case 
(U = 1.595*3;, V = OMx), OTE pairing is indeed the 
most stable as shown in Fig. [SjD). This pairing has s- 
wave gap, whose kx > portion is shown in the inset (k) 
of Fig. EIJD). There are also no nodes of A near kx = kp 
in momentum space, and instead there is a node in Mat- 
subara frequency space as shown in the inset (oj) of Fig. 
[3lJD). It is remarkable that the realization of Berezinskii's 
pairing has been verified based on a microscopic calcula- 
tion. These four cases are summarized in Fig. |4l 

To summarize, we have shown that odd-o; spin-triplet 
s-wave pairing originally proposed by Berezinskii can be 
realized in systems with strong one-dimensionality when 
the strength of charge fluctuation exceeds over that of 
spin fluctuation. Experimentally, it is interesting to look 
for this exotic pairing in QID materials where spin and 
charge fluctuations coexist. Moreover, we may expect 
applying magnetic field to enhance spin-triplet pairing 
(even when charge fluctuation is not so strong) [14;], so 
combining this effect with the strong one-dimensionality 
may increase chances for realizing Berezinskii's pairing 
in actual materials. Also, we hope anomalous properties 
specific to this pairing will be observed in QID systems 
1,16]. 
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FIG. 4: (color online). The most stable pairings dependent 
on the strength of one- dimensionality and spin/charge fluctu- 
ation. 
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